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A  B  BRE  VIATIONS  AND  SYMBOLS 


VK1 
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Rpxq 


plant  matrix 

Input  matrix 
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controller  matrices 

transducer  matrix 
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k-th  Markov  parameter 
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discrete-time  set  {0,T,2T,...} 

set  of  transmission  zeros 

set  of  'slow'  modes 

•«t  of  'fast*  modes 
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1.1  Introduction 


In  spite  of  the  extensive  effort  which  has  been  expended 
during  the  past  twenty  years,  most  currently  available 
techniques  (Athans  and  Falb  1966,  Rosenbrock  1974,  MacFarlane 
1980,  Wonham  1974,  Wolovich  1974,  Davison  and  Ferguson  1981) 
for  the  design  of  tracking  systems  incorporating  linear 
multivariable  plants  are  both  conceptually  and  computationally 
rather  complicated.  Furthermore,  such  techniques  are  almost 
exclusively  concerned  with  the  design  of  analogue  controllers 
and  are  therefore  inapplicable  to  the  much  more  important 
practical  task  of  designing  digital  controllers. 

The  design  techniques  described  in  this  report  are  both 
conceptually  and  computationally  simple.  In  addition,  these 
techniques  are  equally  applicable  to  the  design  of  both 
analogue  and  digital  controllers.  This  universality  and 
simplicity  derives  from  the  fact  that  these  techniques  are 
based  upon  the  systematic  exploitation  of  a  solitary  system- 
theoretic  result  from  the  theory  of  singular  perturbations 
(Porter  and  Shenton  1975) .  This  system- theoretic  result 
exhibits  the  distinctive  asymptotic  structure  of  the  transfer 
function  matrices  of  linear  multivariable  systems  with  'slow* 
and  'fast'  inodes  in  a  manner  which  is  directly  applicable  to 
the  design  of  tracking  systems  incorporating  either  high- 
gain  analogue  controllers  or  fast-sampling  digital  controllers. 

1.2  Regular  and  Irregular  Plants 

In  the  context  of  the  design  of  high-performance  tracking 
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systems,  it  is  convenient  to  classify  linear  multivariable 
plants  as  follows: 

(i)  Regular  plants :  first  Markov  parameter  of  maximal 
rank,  minimum  phase; 

(ii)  Irregular  plants:  first  Markov  parameter  of  non- 
maximal  rank,  non-minimum  phase. 

The  design  of  tracking  systems  incorporating  regular  plants 
is  considered  in  Chapters  2  and  3,  whilst  the  design  of 
tracking  systems  incorporating  irregular  plants  is  considered 
in  Chapters  4  and  5.  In  the  case  of  regular  plants,  arbitrarily 
fast  and  non-interacting  tracking  behaviour  is  achieved: le  by 
appropriate  selection  of  gain  parameters  or  sampling  periods 
simply  by  implementing  error- actuated  analogue  or  digital 
proportional-plus-integral  controllers.  However,  in  the  case 
of  irregular  plants,  good  tracking  behaviour  is  achievable 
only  by  the  simultaneous  implementation  of  inner-loop  com¬ 
pensators  to  provide  appropriate  extra  measurements  for 
control  purposes. 

This  result  is  crucially  important  since  it  makes  explicit 
the  intuitively  obvious  fact  that  controller  and  transducer 
designs  are  inseparable,  and  confirms  the  fact  that  the 
irrelevance  of  much  of  'modern'  control  theory  to  practical 
engineering  design  derives  from  its  failure  to  consider  in¬ 
tegrated  controller/transducer  systems.  The  importance  of 
the  difference  between  regular  and  irregular  plants  obviously 
requires  that  the  set  of  transmission  zeros  of  linear  multi- 
variable  plants  (Portar  and  D'Azzo  1977)  be  readily  computable 
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and  very  powerful  software  for  this  purpose  has  accordingly 
been  developed  (Porter  1979) . 
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CHAPTER  2 


DESIGN  OF  TRACKING  SYSTEMS  INCORPORATING 
HIGH-GAIN  ERROR-ACTUATED  CONTROLLERS 
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2.1  Introduction 


In  this  chapter, singular  perturbation  methods  are  used 
to  exhibit  the  asymptotic  structure  of  the  transfer  function 
matrices  of  continuous -time  tracking  systems  incorporating 
linear  multivariable  plants  which  cure  amenable  to  high-gain 
error- actuated  analogue  control  (Porter  and  Bradshaw  1979a) . 
Such  tracking  systems  consist  of  a  linear  multi variable  plant 
governed  on  the  continuous-time  set  T  ■  [o,+®)  by  state  and 
output  equations  of  the  respective  forms  (Porter  and 
Bradshaw  1979b) 


*x1  ( t) 
x2(t). 


"All 

'  A12 

Xj^tf 

+ 

‘o' 

-A21 

'  A22- 

x2(t). 

-B2- 

and 


y  (t) 


x1(t) 

x2(t) 


f 


(2.1) 


(2.2) 


together  with  a  high-gain  error-actuated  analogue  controller 
governed  on  T  by  a  control- law  equation  of  the  form 


u ( t)  »  g{KQe(t)+K1z(t)  }  (2.3) 

which  is  required  to  generate  the  control  input  vector  u(t) 
so  as  to  cause  the  output  vector  y(t)  to  track  any  constant 
command  input  vector  v(t)  on  T  in  the  sense  that  the  error 
vector  e(t)  ■  v(t)  -  y(t)  assumes  the  steady-state  value 

lim  e(t)  -  lim{v(t)-y(t)}  ■  0  (2.4) 

t-**»  t-*oo 
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for  arbitrary  initial  conditions.  In  equations  (2 . 1) , (2.2) , (2. 3) , 

and  (2.4)  ,xL(  t)  £  lP~l ,  x2(t)ER*,  u(t)  £  R*,  y(t)  £  R*, 

.  (C « (n-£)  x(n-£)  .  £7«(n-i)x£  .  £T  p£x(n-£)  .  _£x£ 

AUCK  /A12^K  ,  A2L^K  ,  A22 K 

B2£RW,  C2eRM,  rank  C2B2  -  t,  eltie^, 

v(t)£R*,  z(t)  =  z(o)  +  £  edt  Cr£(  KqGRAx£, 
and  g  £  R+. 

It  is  evident  from  equations  (2.1)  ,  (2.2)  ,  (2. 3)  that 
such  continuous-time  tracking  systems  are  governed  on  T  by 
state  and  output  equations  of  the  respective  forms 


x^t) 

x2(t) 


9B2Ki  '  A2l“gB2KoCl  '  A22~gB2KoC2^  Lx2 


gfioK 


xL(t) 


(2.5) 


y (t)  -  [o  ,  cx  ,  c2]  Xl(t) 

x2(t). 


(2.6) 


The  transfer  function  matrix  relating  the  plant  output  vector 
to  the  command  input  vector  of  the  closed-loop  continuous¬ 
time  tracking  system  governed  by  equations  (2.5)  and  (2.6)  is 
clearly 
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XIi  *  C1  '  C2 

-1 

'  **  ‘ 

cm-fb.c^cj 

0  '  XIn-i”All  '  “A12 

0 

--gBjKji’Ajj+gBjI^Cj^ .  *I£"A22+9B2KoC2- 

-«B2Ko 

(2,7) 

and  the  high-gain  tracking  characteristics  of  this  sytem  can 
accordingly  be  elucidated  by  invoking  the  results  of  Porter 
and  Shenton  (1975)  from  the  singular  perturbation  analysis 
of  transfer  function  matrices. 

These  results  yield  the  asymptotic  form  of  G(X)  as  the 
gain  parameter  g  -*■  «  and  thus  greatly  facilitate  the  deter¬ 
mination  of  controller  matrices  K  and  K.  such  that  the 

o  l 

continuous -time  tracking  behaviour  of  the  closed-loop  system 
becomes  increasingly  'tight*  and  non-interacting  as  g  is 
increased.  The  frequency-response  and  step-response 
characteristics  of  a  continuous-time  tracking  system  in¬ 
corporating  an  open-loop  unstable  chemical  reactor  (MacFarlane 
and  Kouvaritakis  1977)  are  presented  in  order  to  illustrate 
these  general  results. 

2.2  Analysis 

It  is  evident  from  equation  (2 .7)  thatrby  regarding 
e  *  1/g  as  the  perturbation  parameter,  the  asymptotic  form 
of  the  transfer  function  matrix  GU)  of  the  continuous-time 
tracking  system  as  g  -*•  »  can  be  determined  by  invoking  the 
results  of  Porter  and  Shenton  (1975)  from  the  singular 
perturbation  analysis  of  transfer  function  matrices.  Indeed. 
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these  results  Indicate  that  as  g  ■+■  •  the  transfer  function 
matrix  G{\)  assumes  the  asymptotic  form 


r(A)  *  r (x)  +  f (x ) 


(2.8) 


where 


r  (A) 

r  (A) 


C(  AI-A) _1B 
o  n  o  o 

C2  ( AI^“9A4)  gBjK 


r  , 

LA12C21ko1k1  ' 


A11"A12C21c1- 


B 


LA12C21J 


Co  *  &o\  '  °J 


and 


(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 


A4  "  "B2KoC2  •  (2.14) 

It  follows  from  equations  (2.8)  ,  (2. 9),and(2. 11)  that  the  'slow* 
modes  Zg  of  the  tracking  system  correspond  as  g  +  •  to  the 
poles  Z^^J  Z2  of  ?(A)  where 


•  {A  £C:  |  AK0+K1| «0) 


(2.15) 


and 

z2  *  {1€lciUiI>.t-ji11+»12c;lc1|^»  ,  (2.16) 

and  from  equations (2. 8) , (2. 10) , and (2. 141  that  the  'fast' 
modes  Zf  of  the  tracking  system  correspond  as  g  -  to 


9 


the  poles  Z3  of  r ( X)  where 

Z3  -  U  CCs|XIt+gC2B2Ko|-0}  .  (2.17) 

Furthermore,  It  follows  from  equations  ( 2.9) , (2. 11) , (2. 12) , 
and  (2. 13) that  the  'slow'  transfer  function  matrix 

f ( X)  -  0  (2.18) 

and  from  equation  (2. 10)  and(2.14)  that  the  'fast*  transfer 
function  matrix 

f (X)  *  (XI4+gC2B2Ko)“1gC2B2Ko  .  (2.19) 

Hence,  in  view  of  equations (2.18)  and  (2.19),  it  is 
evident  from  equation  (2.8)  that  asg  +  •  the  trams fer 
function  matrix  G(X)  of  the  continuous  ■‘time  tracking  system 
assumes  the  asymptotic  form 

T(X)  *  (XI4+gC2B2K0)_1gC2B2K0  (2.20) 

in  consonance  with  the  fact  that  only  the  'fast'  modes 
corresponding  to  the  poles  Z3  remain  both  controllable  amd 
observable  as  g  +  ®.  Indeed,  the  'slow*  tramsfer  function 
matrix  T(X)  vanishes  precisely  because  the  'slow'  inodes 
corresponding  to  the  poles  Z^  become  asymptotically  un¬ 
controllable  a*  g  +  »  in  view  of  the  block  structure  of 
the  matrices  AQ  and  BQ  in  equations  (2.11)  and  (2.12)  whilst  the 
'slow'  modes  corresponding  to  the  poles  Z2  become  asymptotically 
unobservable  M  g  +  «  in  view  of  the  block  structure  of  the 
matrices  AQ  and  CQ  in  equations  (2.11)  and  (2.13) 


3.  Synthesis 


It  is  evident  from  equations  (2.5)  and  (2.6)  that  tracking 
will  occur  in  the  sense  of  equation  (4)  provided  only  that 

Z8UZfCc-  (2.21) 

where  C"*  is  the  open  left  half-plane.  In  view  of  equations 
(2.15) , (2.16) ,and(2. 17)£he  'slow*  and  'fast'  modes  will 
satisfy  the  tracking  requirement  (2.21)  for  sufficiently 
large  gains  if  the  controller  matrices  Kq  and  are  chosen 
such  that  both  Z^Cc"  and  Z^Cc"  in  the  case  of  minimum- 
phase  plants  for  which  (Porter  and  D'Azzo  1977)  the  set  of 
transmission  zeros 

Zt  -  (xGc.-AI^^-A^A^^I-OJCC-  (2.22) 

since  it  is  then  immediately  obvious  from  equation  (2.16) 
that  ZjCc".  Moreover,  in  such  cases,  tracking  will  become 
increasingly  'tight'  as  g  -*■<*>  in  view  of  equation  (2.20). 
Furthermore,  if  KQ  is  chosen  such  that 

CjBjKq  m  diag{a^,C2» • . . ,0^}  (2.23) 

where  £1 R+  ( j*l,2, . . . ,1) ,  it  follows  from  equation  (2.20) 
that  the  transfer  function  matrix  G(X)  of  the  continuous¬ 
time  tracking  system  will  assume  the  diagonal  asymptotic 
form 
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and  therefore  that  increasingly  non-interacting  tracking 
behaviour  will  occur  as  g  -►  •  . 


2.4  Illustrative  Example 


These  general  results  can  be  conveniently  illustrated 
by  designing  a  high-gain  error-actuated  analogue  controller 
for  an  open-loop  unstable  chemical  reactor  governed  on  T 
by  the  respective  state  and  output  equations  (MacFarlane 
and  Kouvaritakls  1977} 


xx(t) 

x2(t) 

x3(t) 


1.38  ,  -0.2077  ,  6.715  ,  -5. 

-0.5814  ,  -4.29  ,  0,0. 

1.067  ,  4.273  ,  -6.654  ,  5. 

.  0.048  ,  4.273  ,  1.343  ,  -2.104J  l_x4<t)J 


X^tf 

x2(t) 

x3(t) 

J 

O  ,  0 

5.679  ,  0 

1.136  ,  -3.146 

|_1.136  ,  0 


lu2<t>J 


(2.25) 


and 


Vj.  <t>1 

J1' 

/2(t,J 

* 

o, 

ri ,  o,  i ,  -ii 
[p,i,o,oj 


Xj^t) 

Xj(t) 

x3(t) 

L*4(t). 


(2.26) 


from  which  it  can  be  readily  verified  that  Z.  •  {-1.192, 
-5.039)  Cl  C  and  that  the  first  Markov  parameter 
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C2B2 


'  O  -3.146k 
.  :5.679  .,  O;  J  ; 


(2.27) 


is  of  full  rank.  In  case  -  {1,1}  and  -  2KQ,  it 

follows  from  equations  (2. 3)  ,  (2.23)«and(2.27)that  the  correspond- 
ing  high-gain  error-actuated  analogue  controller  is  governed 
an  by  the  control- law  equation 

-  ’  ■  «■  .  .  ; 

* 

u^t)]  C  0  ,0.1761]  |‘e1(t)l  [  0  ,0.3522]  fz-ft)] 

■g.  + 

u2 <t)J  L**0. 3179 ,  0  JLe2(t)J  [.-0.6358,  0  J  L*2(t)J 

.  .  (2.28) 

and  it  is  evident  from  equations  (2.17)’,  (2.18)  ,  and  (2.19)  that 
Z1  -  {-2,-2} ,■  Z2  *  {-1.192/-5.039},  and  Z 3  -  {-g,-g>.  It 
is  also  evident  from  equation  T2r.  24)  that  the  asymptotic 
transfer  function  matrix  assumes  the  diagonal  form 


'  °1 


o  , 


(2.29) 


and  therefore  that  the  closed-loop  continuous-time  tracking 

'  >  ..  .  S'  »1  ’*  "  J 

system  incorporating  the  chemical  reactor  will  exhibit 

*,  ; '  •.’*'!  v  -  .  •  C-  . .  '  *  • 

increasingly  'tight'  and  non-interacting  tracking  behaviour 
as  g  '♦  -  when  the  "control  input  vector  [u.  (t)  ,  u,(t)]T  is 

*  *>-•  * '  j  •  ■ 

generated  by  the  high-gain  analogue  controller  governed  on 
T  by  equation  (2.28). 

The  actual  frequency-response  loci  G(iu)  for  u  £(-»,+-)  are 
shown  in  Pigs  2.1,  2.2,  and 2. 3  when  g  ■  25,  50,  and  100, 
respectively  and  it  is  clear  that  the  actual  frequency- response 


t  -*■ 


*%»•>* 


loci  approach  tha  asymptotic  frequency-response  loci  F(iu») 
as  tha  gain  paramatar  g  is  incraasad.  Tha  corresponding 
step-response  characteristics  ara  shown  in  Fig  2.4  and  Fig  2.5 
when  [v^t)  ,  Vjlt)]^  ■  [l  ,  0]T  and  [v^(t)  ,  Vjlt)]1*1  ■ 

[o  ,  1] T,  respectively,  and  it  is  avidant  that  incraasingly 
'tight'  and  non- interacting  tracking  occurs  as  tha  gain 
parameter  g  is  increased. 

2.5  Conclusion 

Singular  perturbation  methods  have  bean  used  to  exhibit 
the  asymptotic  structure  of  the  transfer  function  matrices 
of  continuous-time  tracking  systems  incorporating  linear 
multivariable  plants  which  are  amenable  to  high-gain  error- 
actuated  analogue  control.  It  has  bean  shown  that  these 
results  greatly  facilitate  the  determination  of  controller 
matrices  which  ensure  that  tha  closed-loop  behaviour  of  such 
continuous-time  tracking  systems  becomes  increasingly  'tight' 
and  non-interacting  as  the  gain  parameter  g  is  increased. 

These  general  results  have  been  illustrated  by  the  presenta¬ 
tion  of  the  frequency-response  and  step-response  character¬ 
istics  of  a  continuous-time  tracking  system  incorporating 
an  open- loop  unstable  chemical  reactor. 
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Fig  2.4(c) :  g  =  100 
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Fig  2.5(b)  : 


EE 


CHAPTER  3 


DESIGN  OF  TRACKING  SYSTEMS  INCORPORATING 
FAST-SAMPLING  ERROR-ACTUATED  CONTROLLERS 


3.1  Introduction 


In  this  chapter ,  singular  perturbation  methods  are  used 
to  exhibit  the  asymptotic  structure  of  the  transfer  function 
matrices  of  discrete-time  tracking  systems  incorporating 
linear  multivariable  plants  which  are  amenable  to  fast¬ 
sampling  error-actuated  control  (Bradshaw  and  Porter  1980a) . 
Such  tracking  systems  consist  of  a  linear  multivariable 
plant  governed  on  the  continuous -time  set  T  =  [0,+®)  by 
state  and  output  equations  of  the  respective  forms  (Bradshaw 
and  Porter  1980b) 


Xx(t) 

— 

A11  '  A12 

ft 

_ J 

+ 

0 

u(  t) 

_x2  (t) 

-A21  '  A22- 

Lx2(t)J 

-B2- 

(3.1) 


and 


y  (t) 


xL(t) 

x2(t) 


(3.2) 


together  with  a  fast-sampling  error-actuated  digital  con¬ 
troller  governed  on  the  discrete-time  set  TT  =  {0,T,2T,...} 
by  a  control- law  equation  of  the  form 


u(kT)  =  f{Koe(kT)+K1z<kT) }  (3.3) 

which  is  required  to  generate  the  piecewise-constant  control 
input  vector  u(t)  *>  u(kT)  ,  t  £[kT,  (k+1) T)  ,  kT  ^T^.,  so  as 
to  cause  the  output  vector  y(t)  to  track  any  constant  command 
input  vector  v(t)  on  T T  in  the  sense  that  the  error  vector 
e(t)  =*  v(t)  -  y(t)  assumes  the  steady-state  value 
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(3.4) 


lim  e(kT)  =  lim{v(kT) -y (kT) }  ■  0 
k-H»  k-*-<=° 

for  arbitrary  initial  conditions  where  f  =  1/T  and  T  is 

the  sampling  period.  In  equations  (3. 1) , (3. 2) , (3. 3) ,  and  (3.4), 

x1(t)£Rn"*,  x2(t)£R*,  u(t)£R£,  y  (t)  £  R£,  An  £  R(n_£)  x(n"£>  , 

A12£R(n“£)x£,  A21£R£x(n-£>,  a22£r£x£,  B2£r£x£, 

C1£R£x(n_£),  C2£r£x£,  rank  C2B2  =  i,  e(t)£R£,  v(t)  £  R£, 

z (kT)  -  z(o)  +Tj  I  Xe(jT)£R£,  K  £R£x£,  K.  £r£x£,  and 
_  +  j-0  °  1 

f  £r  . 

It  is  evident  from  equations  (3.1) ,  (3.2) ,  and  (3.3)  that 
such  discrete-time  tracking  systems  sure  governed  on  by 
state  and  output  equations  of  the  respective  forms 


z  {(k+1)  T} 
x1(  (k+1)  T) 
x2( (k+1) T) 


and 


■  ' 

-TC1  ,  -TC2 

'  z(kT)  ' 

- 

£*1*1' 

*ll"f'FlK0Cl/*12-f'i'lK0C2 

xx(kT) 

fV2Xl' 

•2rfVoCl'#22-fVoC2- 

x2(kT). 

TI, 


fVo 


v(kT) 


(3.5) 


y  (kT) 


[0  ,  Cj. 


’  z  (kT)  * 
xx(kT) 


(AT) 


(3.6) 
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where 


r 

*11  '  *12 

A11  '  A12 

=  exp' 

T 

,$2i  »  *22- 

k 

-A21  '  A22- 

and 


' 

*1 

fT 

=  J  expl 

-*2- 

O 

A11  '  A12 


,  A 


22-> 


0 

B-» 


dt 


(3.8) 


The  transfer  function  matrix  relating  the  plant  output 
vector  to  the  command  input  vector  of  the  closed-loop 
discrete-time  tracking  system  governed  by  equations  (3.5) 
and  (3.6)  is  clearly 


G(  A) 


'XW  TC1  -  TC2 

-1 

«[0,Cl,C2J 

-M1K1  'XIn-r«ll+£,,lKoCl'-*12«  W2 

ft.  K 

1  o 

--«2K1'  •>Ie-*22+£,,2Koc2. 

f't'-K 
u  2  o-1 

(3.9) 


and  the  fast-sampling  tracking  characteristics  of  this 
system  can  accordingly  be  elucidated  by  invoking  the  results 
of  Porter  and  Shenton  (1975)  from  the  singular  perturbation 
analysis  of  transfer  function  matrices. 

These  results  yield  the  asymptotic  form  of  G(X)  as  the 
sampling  frequency  f  -►  «  and  thus  greatly  facilitate  the 
determination  of  controller  matrices  KQ  and  such  that  the 
discrete-time  tracking  behaviour  of  the  closed-loop  system 
becomes  increasingly  'tight'  and  non-interacting  as  f  is 
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increased.  The  frequency-response  and  step-response 
characteristics  of  a  discrete-time  tracking  system  incorporat¬ 
ing  an  open- loop  unstable  chemical  reactor  (MacFarlane  and 
Kouvaritakis  1977)  are  presented  in  order  to  illustrate  these 
general  results. 


3.2  Analysis 


It  is  evident  from  equation  (3.9)  that,  by  regarding 
e  *  1/f  as  the  perturbation  parameter,  the  asymptotic  form 
of  the  transfer  function  matrix  G(\)  of  the  discrete-time 
tracking  system  as  f  ■*  «  can  be  determined  by  invoking  the 
results  of  Porter  and  Shenton  (1975)  from  the  singular 
perturbation  analysis  of  transfer  function  matrices.  Indeed, 
since  it  follows  from  equations  (3.7)  and  (3.8)  that 


lim  f 

*ll_In-£  '  *12 

A11  '  A12 

f-M9 

H 

t 

Ol 

CM 

O 

r-t 

N 

-A21  '  A22- 

and 


rv 

'o' 

lim  f 

X 

■ 

f-H» 

*2- 

-B2- 

(3.10) 


(3.11) 


these  results  indicate  that  as  f  +  •  the  trams fer  function 
matrix  G(A)  assumes  the  asymptotic  form 

r(A)  -  r (A)  +  r (A)  0.12) 


where 


r<M  -  c0<*VVtV'1tb0 


(3.13) 
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9 


(3.14) 


r  ( X)  **  C2(XIJt-IA-A4)"1B2Kt 


-<*1  '  0 

LA12C2VK1  '  All_A12C2lclJ 


B, 


-A12C21J 


co  *  '  °] 


and 


(3.15) 


(3.16) 

(3.17) 


A.  =■  .  (3.18) 

4  l  O  Z 

It  follows  from  equations  (3.12) , (3. 13) ,  and  (3.15)  that  the 

•slow'  modes  Z  of  the  tracking  system  correspond  as  f  ♦  • 
s 

to  the  poles  Z^U  Z2  of  T(X)  where 

Zx  *  (X  £[C:  j  XI^-I^+TK^KjJ-O}  (3.19) 


and 

Z2  -  (X  €lC:|XIn_A>In_A-TA11+TA12C“1C1|-0}  ,  (3.20) 

and  from  equations  (3. 12) , (3. 14) ,  and  (3.18)  that  the  'fast' 
modes  Zf  of  the  tracking  system  correspond  as  f  -►  •  to  the 

A 

poles  Z3  of  r ( X)  where 

Z3  -  U  £.C: | XIa-Ia+C2B2K0 |-0}  .  (3.21) 

Furthermore,  it  follows  from  equations  (3.13) , (3.15) , (3.16) , 
and(3.17)  that  the  'slow'  transfer  function  matrix 


r(X)  -  0 


(3.22) 
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and  from  equations  (3.14)  and  (3.18)  that  the  'fast1  transfer 
function  matrix 

f(X)  -  (XiJt-I1+C2B2K0)“1C2B2K0  .  (3.23) 

Hence,  in  view  of  equations  (3.22)  and  (3.23),  it  is  evident 
from  equation  (3. 12)  that  as  f  -*•  «  .  the  transfer  function 
matrix  G(X)  of  the  discrete-time  tracking  system  assumes 
the  asymptotic  form 

r  ( X)  -  (Xl4-IA+C2B2Kor1C2B2K0  ~  (3.24) 

in  consonance  with  the  fact  that  only  the  ’fast1  modes 
corresponding  to  the  poles  remain  both  controllable  and 
observable  as  f  +  •.  Indeed,  the  'slow*  transfer  function 
matrix  T(X)  varnishes  precisely  because  the  'slow1  modes 
corresponding  to  the  poles  Z1  became  asymptotically  un¬ 
controllable  as  f  •  in  view  of  the  block  structure  of  the 
matrices  AQ  and  BQ  in  equations  (3.15)  and  (3.16)  whilst  the 
'slow'  modes  corresponding  to  the  poles  Z2  become  asymptotically 
unobservable  as  f  -*■  »  in  view  of  the  block  structure  of  the 
matrices  kQ  and  CQ  in  equations  (3.15)  and  (3.. 7). 

3.3  Synthesis 

It  is  evident  from  equations  (3.5)  and  (3.6)  that  tracking 
will  occur  in  the  sense  of  equation  (3.4)  provided  only  that 

ZBUZfO"  (3.25) 

where  0~  is  the  open  unit  disc.  In  view. of  equations  (3.19) 
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(3.20),and(3.21)  ,the  'slow*  and  'fast'  modes  will  satisfy 
the  tracking  requirement  (3.25)  for  sufficiently  fast  sampling 
frequencies  if  the  controller  matrices  Kq  and  are  chosen 
such  that  both  Z^d.'D  for  sufficiently  small  sampling 
periods  and  in  the  case  of  minimum-phase  plants  for 

which  (Porter  and  D'Azzo  1977}  the  set  of  transmission  zeros 

Zfc  -  U  CCsIXIh^-Ajj+A^C^CjJ-O}  C.C"  (3.26) 

where  C~  is  the  open  left  half-plane  since  it  is  then 

immediately  obvious  from  equation  (3. 20) that  Z2dv~  for 

sufficiently  small  sampling  periods.  Moreover ,  in  such  cases, 

tracking  will  become  increasingly  'tight1  as  f  -►  <*>  in  view 

of  equation  (3.24) .  Furthermore,  if  K  is  chosen  such  that 

o 

^"2® 2^o  **  ,(*2 '  *  *  *  i,  ^  .(3.27) 

where  l-o^  £  R  H  P  (j«l,2, ...,£) ,  it  follows  from  equation 
(3.24) that  the  transfer  function  matrix  G(X)  of  the  discrete-time 
tracking  system  will  assume  the  diagonal  asymptotic  form 

r(X)  -  diag{x_1+ai  »  x-l+a2  •  •••  ’  X-l+cr^}  (3*28) 

and  therefore  that  increasingly  non-interacting  discrete¬ 
time  tracking  behaviour  will  occur  as  f  -►  •. 

3.^  Illustrative  Example 

These  general  results  can  be  conveniently  illustrated 
by  designing  a  fast-sampling  error-actuated  digital  controller 
for  an  open-loop  unstable  chemical  reactor  governed  on  T 
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by  the  respective  state  and  output  equations  (MacFarlane 
and  Kouvaritakis  1977) 


"x^t) 

x2(t) 

X3(t) 

x4(t)_ 

L 

1.38  ,  - 0.20 1 

-0.5814  ,  -4.29 
1.067  ,  4.27: 

0.048  ,  4.27: 

0  ,  o' 

5.679  ,  0 

1.136  ,  -3.146 

1.136  ,  0 


u^t) 

«2(t)- 


Xj^  ( t) 

x2(t) 

X3(t) 

u 

-X4(t). 

(3.29) 


and 


y  ( t) 


1  ,  0  ,  1  ,  -1 
i_0  ,  1  ,  0  ,  0. 


Xx(t) 

*2<« 

X3(t) 

x4(t)J 


(3.  30) 


from  which  it  can  be  readily  verified  that  Zfc  a  (-1.192, 
-5.039)  Cl  C  and  that  the  first  Markov  parameter 


C2®2 


0 

.5.679 


-3.146 


(3.  31) 


In  case  {a^#cr2}  “  {1/1}  end  •  2XQ,  it  follows  from 
equations  (3.3) ,  (3. 2 7), and (2. 38) that  the  corresponding  fast 
sampling- error-actuated  digital  controller  is  governed  on 
Tt  by  the  control-law  equation 
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lr 

'ux(kT) 

-f 

u2(kT) 

,0.1761 

,  0 


e^(kT) 

e2(kT) 


0  ,0.3522 

0.6358,  0 


*x(kT) 

z2  (kT) 


(3.32) 


and  It  is  evident  from  equations  (3. 19) , (3.20) ,  and  (3.21)  that 
«  U-2T,1-2T>,  Z2  -  (l-l. 192T, 1-5.039T) ,  and  Z3  =  (0,0). 

It  is  also  evident  from  equation  (3.24)  that  the  asymptotic 
transfer  function  matrix  assumes  the  diagonal  form 


r  (X) 


(3.33) 


and  therefore  that  the  closed-loop  discrete-time  tracking 
system  incorporating  the  chemical  reactor  will  exhibit 
increasingly  'tight'  and  non-interacting  tracking  behaviour 
as  f  -*■  •  when  the  piecewise-constant  control  input  vector 
[ux(t)  ,  u2(t)]T  -  [u^(kT)  ,  u2(kT)]T,  t£[kT  ,  (k+l)T)  , 
kT  fT,  is  generated  by  the  fast-sampling  digital  con¬ 
troller  governed  on  by  equation  (3.32). 

The  actual  frequency-response  loci  G(eiwT)  for  o>T  ,  2ir] 

are  shown  in  Figs  3.1,  3.2,  and  3.3  when  T  *0.04  ,  0.02,  and  0.01, 
respectively,  and  it  is  clear  that  the  actual  frequency- 
response  loci  approach  the  asymptotic  frequency-response 
loci  r(eluT)  as  the  sampling  frequency  f  is  increased.  The 
corresponding  step-response  characterietics  are  shown  in 
Pig  3. 4  and  Fig  3.5  when£v^(t)  ,  V2(t>]T  -  [l  ,  o]T  and  [v^t)  , 
v2(t)]  »  [o  ,  l]T,  respectively,  and  it  is  evident  that  in¬ 
creasingly  'tight'  and  non-interacting  tracking  occurs  as  the 
sampling  frequency  f  is  increased. 
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3,5  Conclusion 


Singular  perturbation  methods  have  been  used  to  exhibit 
the  asymptotic  structure  of  the  transfer  function  matrices 
of  discrete-time  tracking  systems  incorporating  linear  multi- 
variable  plants  which  are  amenable  to  fast-sampling  error- 
actuated  digital  control.  It  has  been  shown  that  these 
results  greatly  facilitate  the  determination  of  controller 
matrices  which  ensure  that  the  closed- loop  behaviour  of 
such  discrete-time  tracking  systems  become  increasingly 
'tight'  and  non-interacting  ae  the  sampling  frequency  f  is 
increased.  These  general  results  have  been  illustrated  by 
the  presentation  of  the  frequency-response  and  step-response 
characteristics  of  a  discrete-time  tracking  system  incorporat¬ 
ing  an  open- loop  unstable  chemical  reactor. 
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0.02 


0.01 


Fig  3.5(b)  : 


Fig  3.5(c) :  T  =  0.01 


CHAPTER  A 


DESIGN  OF  TRACKING  SYSTEMS 
INCORPORATING  INNER-LOOP  COMPENSATORS 
AND  HIGH-GAIN  ERROR-ACTUATED  CONTROLLERS 
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4.1  Introduction 


In  this  chapter,  singular  perturbation  methods  are  used 
to  exhibit  the  asymptotic  structure  of  the  transfer  function 
matrices  of  continuous-time  tracking  systems  incorporating 
linear  multivariable  plants  which  are  amenable  to  high-gain 
error-actuated  analogue  control  (Porter  and  Bradshaw  1979) 
only  if  extra  plant  output  measurements  sure  generated  by 
the  introduction  of  appropriate  transducers  and  processed 
by  inner-loop  compensators.  Such  tracking  systems  consist 
ttif'a  linear  multi variable  plant  governed  on  the  continuous¬ 
time  set  T  **  [o,+»)  by  state,  output,  and  measurement  equations 
of  the  respective  forms 


*x(t)- 

B 

A11  '  A12 

’x^tf 

+ 

r°’ 

u(t) 

*2(t)- 

-A21  '  A22- 

-x2(t). 

LB2- 

y(t) 


Xj^t) 

x2(t)_ 


and 


(4.1) 


(4.2) 


w(t) 


Xj.ttf 

x2(t)_ 


(4.3) 


together  with  a  high-gain  error-actuated  analogue  controller 
governed  on  T  by  a  control-law  equation  of  the  form 


u(t)  -  g{KQe(t)+K1z(t)}  (4.4) 

which  is  required  to  generate  the  control  input  vector  u(t) 
so  as  to  cause  the  output  vector  y(t)  to  track  any  constant 
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command  input  vector  v(t)  on  T  in  the  sense  that 


lim{v(t) -y (t) }  *  0  (4.5) 

as  a  consequence  of  the  fact  that  the  error  vector  e(t)  = 
v(t)  -  w(t)  assumes  the  steady-state  value 

lim  e(t)  =  lim{v(t) -w(t)  }  *  0  (4.6) 

t-*~00  t-H*> 

for  arbitrary  initial  conditions,  tn  equations  (4.1),  (4.2), 

(4.3)  ,and(4.4)^c1(t)  61  x2(t)  61  R£ ,  u(t)  61  R£,  y(t)  ££  R£, 

w(t)£^  A116R(n‘l,x(n’1),  A126j(Hlxi,  a216r1xM), 
a22grm,  b26rm,  ClGRix(n"£),  c2  eR£x£,  FleR£x(n"£) , 

F2  6Ir£x£,  rank  C2b2  <  £,  rank  F2B2  *  £,  e(t)  6R1;  v(t)  61  R£, 
z(t)  -  z (o)  *  /*  edt  61  R£ t  Kq£rM,  Ki£r£x£,  geR+,  and 

[Fj_  f  —  £c^+MAj^  ,  C2+MA^2^]  (4.7) 

where  M  £~R£x^n  ^  •  It  is  evident  from  equations  (4.2) ,(4.3) ,  and 
(4.7)  that  the  vector 


w(t)  -  y(t) 


*l(t) 

x2(t) 


(4.8) 


of  extra  measurements  is  such  that  v(t)  and  y(t)  satisfy 
the  tracking  condition (4. 5)  for  any  MGlR£x^n”£*  if  e(t) 
satisfies  the  steady-state  condition  (4.6)  since  equation  (4.1) 
clearly  implies  that 


lim[Au 


A12l 


x1(t)‘ 

.x2(t). 


(4.9) 
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in  any  steady  state.  However,  the  condition  that  rank  f2B2«  l 
requires  that  C2  £  RAx*  and  A12  are  such  that 

M  can  be  chosen  so  that 

rank  F2  ■  rank(C2+MA12)  ■  A  .  (4.10) 

It  is  evident  from  equations  (4. 1) , (4.2) ,  (4. 3) ,  and  (4.4) 
that  such  continuous- time  tracking  systems  are  governed  on 
T  by  state  and  output  equations  of  the  respective  forms 


’ z(t) " 

y (t)  -  [0  ,  c1  ,  c2]  xL(t)  .  (4.  12) 

-*2(t)- 

The  transfer  function  matrix  relating  the  plant  output 
vector  to  the  command  input  vector  of  the  closed-loop 
continuous-time  tracking  system  governed  by  equations  (4.11) 
and  (4.12)  is  clearly 


G<i)-(b,Cl,C2]  0  '  XI„.t-Au  *  ~k12 

-9B2*1'-A21+982Vl'JIt-A22^B2*op 
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(4-13) 


and  the  high-gain  tracking  characteristics  of  this  system 
can  accordingly  be  elucidated  by  invoking  the  results  of 
Porter  and  Shenton  (1975)  from  the  singular  perturbation 
analysis  of  transfer  function  matrices. 

These  results  yield  the  asymptotic  form  of  G(X)  as 

the  gain  parameter  g  ■*  °°  and  thus  greatly  facilitate  the 

determination  of  controller  and  transducer  matrices  K  ,  K, , 

o  l 

and  M  such  that  the  continuous-time  tracking  behaviour  of 
the  closed-loop  system  becomes  increasingly  non-interacting 
as  g  is  increased.  The  frequency-response  and  step-response 
characteristics  of  a  continuous- time  flight-control  system 
for  the  longitudinal  dynamics  of  an  aircraft  (Kouvaritakis, 
Murray,  and  MacFarlane  1979)  are  presented  in  order  to 
illustrate  these  general  results. 

4.2  Analysis 

It  is  evident  from  equation (4. 13) that,  by  regarding 
e  =  1/g  as  the  perturbation  parameter,  the  asymptotic  form 
of  the  transfer  function  matrix  G(X)  of  the  continuous-time 
tracking  system  as  g  +  •  can  be  determined  by  invoking  the 
results  of  Porter  and  Shenton  (1975)  from  the  singular 
perturbation  analysis  of  transfer  function  matrices.  Indeed, 
these  results  indicate  that  as  g  +  •  the  transfer  function 
matrix  G(X)  assumes  the  asymptotic  form 

r(X)  -  f (X)  +  f (X)  (4.  14) 


where 
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r  (X) 


9 


(4.15) 


and 

Z2  -  (X  £C:  IXI^-Aj^+A^F^fJ-O)  ,  (4.22) 

and  from  aquations  (4.14),  (4.16),  and  (4.20)  that  the  'fast' 
modes  Zj  of  the  tracking  system  correspond  as  g  *•>  *  to  the 

a 

poles  Z3  of  r ( X)  where 

Z3  -  (X  CCslXI^PjBj^l-O)  .  (4.  23) 

Furthermore,  It  follows  from  equations  (4.15),  (4.17),  (4.18), 
and(4. 19)  that  the  'slow'  transfer  function  matrix 


and  from  equations (4. 16)  and( 4. 20) that  the  'fast'  transfer 
function  matrix 


r(X)  =  C2F21(XIJl+gF2B2K0)_1gF2B2K0  .  (4.25) 

Hence,  in  view  of  equations  (4.24)  and  (4.25),  it  is 
evident  from  equation(4. 14)  that  as  g  +  «  the  transfer 
function  matrix  G <  X)  of  the  continuous-time  tracking  system 
assumes  the  asymptotic  form 


T(A)  -  (C1-C2F2  F1) ( XIn_ff-A1 1+A1 2F2  F.,)  Ai2F21 


1  2  2  ‘1'  n-£  11  12  2  *1J 

2B2Ko 


+  C2F21(XI£+gF2B2KQ)'1gF,B,K 


(4.26) 


in  consonance  with  the  fact  that  both  the  ’slow'  modes 
corresponding  to  the  poles  Z2  and  the  'fast'  modes  correspond¬ 
ing  to  the  poles  Z3  possibly  remain  both  controllable  and 
observable  as  g  +  •.  However,  the  'slow'  transfer  function 
matrix  T(X)  reduces  to  the  form  expressed  by  equation  (4.24) 
precisely  because  the  'slow'  modes  corresponding  to  the  poles 
definitely  become  asymptotically  uncontrollable  as  g  -*■  ® 

in  view  of  the  block  structure  of  the  matrices  A  and  B  in 

o  o 

equations  (4.17)  and  (4.18). 


4.3  Synthesis 


It  is  evident  from  equations  (4.7) , (4.9) , (4.11) ,  and  (4.12) 
that  tracking  will  occur  in  the  sense  of  equation  (4.5) 
provided  only  that 


z,Uz(0' 


(4.27) 


where  C~  is  the  open  left  half-plane.  In  view  of  equations (4. 21) , 
(4.22),  and  (4.23),  the  'slow1  and  'fast*  modes  will 
satisfy  the  tracking  requirement (4. 27) for  sufficiently  large 
gains  if  the  controller  and  transducer  matrices  KQ,  K^,  and 
M  are  chosen  such  that  Z^CZc”,  Z2  (Zc”,  and  Z3CZc”  in 
the  case  of  plants  for  which  the  transducer  matrix  m  can  be 
simultaneously  chosen  so  as  to  satisfy  the  measurement  con¬ 
dition  expressed  by  equation  (4.10).  Moreover,  if  Kq  and  M 

A 

are  chosen  so  that  both  r(X)  and  r(X)  are  diagonal  transfer 
function  matrices  by  requiring  that 

F2B2Ko  *  (C2+MA12)B2Kc  -  diag{o1,a2,...,oJl)  (4.28) 

, 

where  Oj£lR  (j»l,2, . . .  ,1)  in  the  case  of  r(X),  it  follows 
from  equation (4. 20) that  the  transfer  function  matrix  G(X) 
of  the  continuous-time  tracking  system  will  assume  the 
diagonal  asymptotic  form  T(X)  and  therefore  that  increasingly 
non-interacting  tracking  behaviour  will  occur  as  g  ♦  «. 
Furthermore,  such  tracking  behaviour  will  exhibit  high 
accuracy  in  the  face  of  plant-parameter  variations  provided 
that  the  steady-state  conditions  expressed  by  equation  (4.9) 
correspond  to  'kinematic'  relationships  which  hold  between 
the  state  variables  as  a  consequence  of  the  fundamental 
dynamical  structure  of  the  plant. 

However,  increasingly  'tight*  tracking  behaviour  will 
not  in  general  occur  as  g  ■*>  •  in  view  of  the  possible  presence 
in  the  plant  output  vector  of  'alow'  modes  corresponding  to 
the  poles  Z2  of  the  'alow'  transfer  function  matrix  T(X). 
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But  the  presence  of  any  such  'slow*  modes  is  the  inevitable 
consequence  of  the  introduction  of  appropriate  transducers 
which  generate  extra  plant  output  measurements  so  as  to 


ensure  that  rank  F^B,.  =  £  and  thus  render  plants  for  which 
rank  C2B2  <  *  a™®114113!®  to  high-gain  error-actuated  control. 
Indeed,  increasingly  'tight'  tracking  behaviour  is  achievable 
as  g  -*•  ®  only  in  the  case  of  plants  for  which  rank  C2b2  -  £ 

(Porter  and  Bradshaw  1979}  and  which  are  accordingly  amenable 
to  high-gain  error-actuated  control  without  the  necessity 
for  the  generation  of  extra  plant  output  measurements. 

4,4  Illustrative  Example 


These  general  results  can  be  conveniently  illustrated 
by  designing  a  high-gain  error-actuated  analogue  flight 
controller  for  the  longitudinal  dynamics  of  an  aircraft 
governed  on  T  by  the  respective  state  and  output  equations 
(Kouvaritakis,  Murray,  and  MacFarlane  1979) 
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Xj^t) 

x2  (t) 

*3<t) 

x4(t) 

ml 

Xglt^ 

0.1575  ,  0  ,  -0.0732 


u^t) 

u2(t) 

Lu3  (t)J 


(4.29) 
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and 


1/0,  0,0,0 

y2(t) 

- 

0,0,  1,0,0 

^3(t). 

0  ,  1,0,  0,0^ 

"x1  (t)" 
x2  (t) 
X3  (t) 
x4(t) 
x5(t). 


(4.30) 


from  which  it  can  be  readily  verified  that  the  first  Markov 
parameter 


C2B2 


0  ,  0  ,  o' 

-0.0012  ,1,0 

O  ,0,0 


(4.  31) 


is  rank  defective.  In  case  {o^o^o-j}  *  {1,1,1},  Kj  *  Kq, 
and 


0.25  ,  O  I 


M 


0  ,  0 
0  ,  0.25 


(4.  32) 


it  follows  from  equations  (4.3),  (4.4),  and  (4.28)  that  the 
corresponding  transducers  and  high-gain  error- actuated 
analogue  controller  are  governed  on  T  by  the  respective 
measurement  and  control- law  equations 


Wj^  ( t)' 

"l  ,  0.283 

f  0 

,  0 

,  0.25" 

m 

0  , 

0 

,  1 

,  0 

,  0 

.w3(t) 

p  , 
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,  0 

,  0.25 

/  0  . 

x^t  )* 
x2(t) 
x3(t) 
x4(t) 

„x5(t)- 


(4.33) 
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and 


' 

v  — T 

e 

(-* 

ft 

-28.97 

,  0  ,  -1.274 

e1(t> 

u2(t) 

*  9‘ 

-0.0348 

,  1  ,  -0.0015 

e2(t) 

u3(t). 

_  -7.690 

,  0  ,  -2.740 

e3(t). 

'-28.97 

,  0  ,  -1.274  ‘ 

zL(t) 

+ 

-0.0348 

,  1  ,  -0.0015 

z2  (t) 

_  -7.690 

,  0  ,  -2.740  _ 

.z3(t)_ 

and  it  is  evident  from  equations  (4.21),  (4.22),  and  (4.23) 
that  -  {-1,-1, -1} ,  Z2  =  {-4,-4},  and  =  {-g,-g,-g). 

It  is  also  evident  from  equation (4. 26) that  the  asymptotic 
transfer  function  matrix  assumes  the  diagonal  form 


ru)  » 


(4.35) 


and  therefore  that  the  closed-loop  continuous-time  flight- 
control  system  for  the  longitudinal  dynamics  of  the  aircraft 
will  exhibit  increasingly  non-interacting  tracking  behaviour 
a*  g  ♦  ■  when  the  control  input  vector  £u^(t)  ,  u2(t)  ,  u3(t)] 
is  generated  by  the  high-gain  analogue  controller  governed 
on  T  by  equation  (4. 34).  However,  it  is  apparent  from  equation 
(4.35)that  increasingly  'tight'  tracking  behaviour  will  be 
achieved  only  in  the  case  of  y2(t)  in  view  of  the  presence 
in  y^(t)  and  y3(t)  of  '.he  'slow1  modes  corresponding  to  the 
poles  2j. 
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The  actual  frequency-response  loci  G(iu>)  for  w  ^ (-*/+«) 
are  shown  in  Figs  4 . 1 , 4 . 2 ,  and  4 . 3  when  g  ■  10,20 ,  and  50 , 
respectively,  and  it  is  clear  that  the  actual  frequency- 
response  loci  approach  the  asymptotic  frequency-response 
loci  r(iu)  as  the  gain  parameter  g  is  increased.  The 
corresponding  step-response  characteristics  are  shown  in 
Fig  4.4,Fig4.5,andFig4.6  vrtien[v1(t)  ,  v2(t)  ,  v3(t)]T  * 
[1,0,  0]T,  [vx(t)  ,  v2(t)  ,  v3(t)]T  »  [0  ,  1  ,  0]T,  and 
[v^ft)  ,  v2(t)  ,  v3(t)]T  “  [o  ,  0  ,  l]T,  respectively,  and 
it  is  evident  that  increasingly  non-interacting  tracking 
occurs  as  the  gain  parameter  g  is  increased  but  that  in¬ 
creasingly  'tight'  tracking  occurs  only  in  the  case  of  y2(t). 

4.5  Conclusion 

Singular  perturbation  methods  have  been  used  to  exhibit 
the  asymptotic  structure  of  the  transfer  function  matrices 
of  continuous-time  tracking  systems  incorporating  linear 
multivariable  plants  which  are  amenable  to  high-gain  error- 
actuated  analogue  control  only  if  extra  plant  output  measure¬ 
ments  are  generated  by  the  introduction  of  appropriate  trans¬ 
ducers  and  processed  by  inner- loop  compensators .  It  has  been 
shown  that  these  results  greatly  facilitate  the  determination 
of  controller  and  transducer  matrices  which  ensure  that  the 
closed-loop  behaviour  of  such  continuous-time  tracking  systems 
becomes  increasingly  non-interacting  as  the  gain  parameter  g 
is  increased.  These  general  results  have  been  illustrated  by 
the  presentation  of  the  frequency-response  and  step-response 
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characteristics  of  a  continuous -time  flight-control  system 
for  the  longitudinal  dynamics  of  an  aircraft. 
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Fig  4.2(e):  G^dtu)  ;  g  =  20 
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Fig  4.4(c): 


Fig  4.5(b)  : 


Fig  4.5(c)  : 


Fig  4. 6(a) 


Fiq  4, 6(b) : 


Pi  q  4.6(c): 


CHAPTER  5 


DESIGN  OF  TRACKING  SYSTEMS 
INCORPORATING  INNER-LOOP  COMPENSATORS 
AND  FAST-SAMPLING  ERROR-ACTUATED  CONTROLLERS 
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5,1  Introduction 


In  this  chapter  singular  perturbation  methods  are  used 
to  exhibit  the  asymptotic  structure  of  the  transfer  function 
matrices  of  discrete-time  tracking  systems  incorporating 
linear  multivariable  plants  which  are  amenable  to  fast- 
sampling  error-actuated  digital  control  (Bradshaw  and  Porter 
1980)  only  if  extra  plant  output  measurements  are  generated 
by  the  introduction  of  appropriate  transducers  and  processed 
by  inner-loop  compensators.  Such  tracking  systems  consist 
of  a  linear  multivariable  plant  governed  on  the  continuous- 
time  set  T  =  [o ,  +°°)  by  state,  output,  and  measurement 
equations  of  the  respective  forms 


and 


"x  i  ( t ) 

x2(t)J 


A11  '  A12 
•A21  '  A22- 


x1(t) 
•x2  * 


lB2J 


u  ( t ) 


(5.1) 


y (t)  «  [cx  ,  c2] 


xi(t) 

Lx2(t)- 


(5.2) 


w(t)  =  [FL  ,  F2] 


xx(t) 
_x2  ( t) 


(5.3) 


together  with  a  fast-sampling  error-actuated  digital  con¬ 
troller  governed  on  the  discrete-time  set  TT  =  {o,T,2T,...} 
by  a  control-law  equation  of  the  form 


u(kT)  -  f{Koe(kT)+K1z(kT)  }  (5.4) 

which  is  required  to  generate  the  control  input  vector 
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so  as  to  cause  the 


u ( t)  *  u(kT)  /  t  CCkT  '  (k+DT),  kT  £Tt, 
output  vector  y(t)  to  track  any  constant  command  input 
vector  v(t)  on  T T  in  the  sense  that 


lim{v(kT)-y(kT) }  «  0  (5.5) 

k-t-oo 

as  a  consequence  of  the  fact  that  the  error  vector  e(t)  = 
v ( t)  -  w(t)  assumes  the  steady-state  value 

lim  e(kT)  »  lim{ v(kT) -w(kT) >  *  0  (5.6) 

k-*»  k-***> 


for  arbitrary  initial  conditions.  In  equations  (5.1),  (5.2), 
(5. 3)pnd  (5. 4)  ,xx  ( t)  £”  Rn_£ ,  x2(t)^RA,  u(t)££R£,  y(t)  (£  R£, 

w(t)eR\  a11£rM)xM)(  A19GR(n“a)xi,  a,,  eR£x(n_i), 


12 


21 


a22£)!m,  b2£r*x\  Cl£R£x(n“l),  c2£rw,  Fi  e  R*x(n_£) , 

F2  £1rAxA,  rank  C_B2  <  i,  rank  F2B2  =•  l,  e(t)££R*,  v(t)  ££  R^, 

z(kT)  »  z(o)  +  Tj  r  e<jT)£lR£,  k Q£Rlxi,  k  £  rAx£,  fD+, 

j-0  o  i. 

and 


[fx  ,  F2]  *  [Cj+MAj^  ,  C2+MA12]  (5.7) 

where  M£2RAx^n”^.  It  is  evident  from  equations  (5.2),  (5.3), 
and  (5.7)  that  the  vector 


w(t)  -  y (t) 


xx(t) 

x2(t) 


(5.8) 


of  extra  measurements  is  such  that  v(kT)  and  y(kT)  satisfy 
the  tracking  condition (5. 5)  for  any  M  £1  jj^xtn-i)  ^  e(kT) 
satisfies  the  steady-state  condition  (5.6)  since  equation  (5.1) 
clearly  implies  that 
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lim[A 

t+oo 


11  '  12 


'x,  (t) 


*2(t) 


(5.9) 


in  any  steady  state.  However,  the  condition  that  rank 


F2B2  =  i  requires  that  C2£LR*'xJ'  and  A^2 
such  that  M  ^  can  ’De  chosen  so  that 


(n-l) xl 


rank  F2  =  rank(C2+MA^2)  =  l 


(5.10) 


It  is  evident  from  equations  (5. 1)  ,  (5.2) , (5. 3)  ,  and  (5.4) 
that  such  discrete-time  tracking  systems  are  governed  on  T 
by  state  and  output  equations  of  the  respective  forms 


"  z {  (k+1)  T}  * 
xL{  (k+1)  T} 
_x2(  (k+1) T }_ 


“  f'i,lKl,<I,ll_f7lKoFl'$12“f'i,lKoF2  xl(kT) 

.fH,2Ki,<S2i-fVoFl'<&22-f'1'2KoF2JLx2(kT) 


z(kT) 


+  fV^K  v(kT) 

-fVo- 


(5.11) 


y(kT)  =  [0  ,  C1  ,  C2]  xx(kT) 

_x2(kT) 


z  (kT) 


(5.12) 


where 
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1 


(5.13) 


' 

r 

$11  '  $12 

A11  '  A12 

=  exp- 

T 

_$2i  '  $22. 

. 

-A2 1  '  A22- 

and 


) 

,T 

An  '  A12 

A 

'o' 

=  J  exp- 

-2- 

V 

> 

A2i  #  A22_ 

J 

-B2- 

(5.14) 


The  transfer  function  matrix  relating  the  plant  output 
vector  to  the  command  input  vector  of  the  closed- loop 
discrete-time  tracking  system  governed  by  equations  (5.11) 
and  (5.12)  is  clearly 


G(  A)  = 


'"r1.-  TPi  -  tp2  i 

-1 

'TIi ' 

Zo,cvc2] 

->i2+£',iKoF2 

fVo 

.-n&v  •XIrt22tfT2KoF2- 

-fVo 

(5.15) 

and  the  fast-sampling  tracking  characteristics  of  this  system 
can  accordingly  be  elucidated  by  invoking  the  results  of 
Porter  and  Shenton  (1975)  from  the  singular  perturbation 
analysis  of  transfer  function  matrices. 

These  results  yield  the  asymptotic  form  of  G ( X )  as  the 
sampling  frequency  f  -*■  ®  and  thus  greatly  facilitate  the 
determination  of  controller  and  transducer  matrices  KQ,  k1# 
and  M  such  that  the  discrete-time  tracking  behaviour  of  the 
closed-loop  system  becomes  increasingly  non-interacting  as 
f  is  increased.  The  frequency-response  and  step-response 
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characteristics  of  a  discrete-time  flight-control  system 
for  the  longitudinal  dynamics  of  an  aircraft  (Kouvaritakis , 
Murray,  and  MacFarlane  1979)  are  presented  in  order  to 
illustrate  these  general  results. 


5.2  Analysis 


It  is  evident  from  equation  (5.15)  that,  by  regarding 
e  =  1/f  as  the  perturbation  parameter,  the  asymptotic  form 
of  the  transfer  function  matrix  G(X)  of  the  discrete-time 
tracking  system  as  f  +  •  can  be  determined  by  invoking  the 
results  of  Porter  and  Shenton  (19  7  5)  from  the  mgular 
perturbation  analysis  of  trains fer  function  matrices.  Indeed, 
since  it  follows  from  equations  (5.13)  and  (5.14)  that 


lim  f 

*ll-In-«.  '  *12 

_ 

A11  '  A12 

f -.CO 

-  *21  '  *22-I£- 

-A21  '  A22- 

and 


lim  f 


<5.16) 


(5.17) 


these  results  indicate  that  as  f  +  «  the  transfer  function 
matrix  G(A)  assumes  the  asymptotic  form 

r( a)  =  r(A)  +  ?(x)  (5.18) 


where 

r ( A)  =  C  ( AI  -I  -TA  ’ -1TB 
o  n  n  o  o 


(5.19) 
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9 


(5.20) 


(5.21) 

(5.22) 

(5.23) 


-  -B2KoF2 


(5.24) 


It  follows  from  equations  (5.18)/  (5.19)/  and  (5.21)  that  the 

'slow'  modes  Z  of  the  tracking  system  correspond  as  f  +  “ 
s 

to  the  poles  Z ^  of  T (X)  where 

Zl  »  (A  €lC:  |  xI4-Ijl+TK“1K1l=0}  (5.25) 


and 


Z2  -  U  €IC:  IAI^-I^-TAj^-KTA^F^fJ-O}  /  (5.26) 

and  from  equations  (5.8)/  (5.20)/  and  (5.24)that  the  'fast' 
modes  Zf  of  the  tracking  system  correspond  as  f  +  ®  to 

A 

the  poles  Z3  of  r ( X )  where 

Z3  »  (X  G  C:  I  XIJl-'Ijl+F2B2Ko|=0}  .  (5.27) 

Furthermore/  it  follows  from  equations  (5.19) ,  (5.21) ,  (5. 22) , 
and (5. 23) that  the  'slow1  trams fer  function  matrix 


(C1-C2F2  fi)  (XIn-JrIn-jTTA114'rA 


F~^ 
12*  2 
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(5.28) 


and  from  equations (5. 20) and (5. 24) that  the  ’fast'  transfer 
function  matrix 


r  (X) 


C2F21(XVVF2B2V"lF2B2Ko 


(5.29) 


Hence,  in  view  of  equations  (5.28)  and  (5.29),  it  is 
evident  from  equation (5. 18) that  as  f  +  »  the  transfer 
function  matrix  G(X)  of  the  discrete-time  tracking  system 
assumes  the  asymptotic  form 


r(X)  =  UIn.1-In.rTA11+M12F-1F1)-1TS12Fj1 


+  C-F“1(XI0-I0+F0B0K  ) -1F0BoK„ 
22  ££  2  2  0  2  2  O 


(5.  30) 


in  consonance  with  the  fact  that  both  the  'slow'  modes 
corresponding  to  the  poles  and  the  'fast'  modes  correspond¬ 
ing  to  the  poles  Z3  possibly  remain  both  controllable  and 
observable  as  f  +  «.  However,  the  'slow'  transfer  function 
matrix  T(X)  reduces  to  the  form  expressed  by  equation  (5.28) 
precisely  because  the  'slow'  modes  corresponding  to  the 
poles  definitely  become  asymptotically  uncontrollable  as 

f  -*■  <»  in  view  of  the  block  structure  of  the  matrices  A  and 

o 

Bq  in  equations  (5.21)  and  (5.22)  . 


5.3  Synthesis 

It  is  evident  from  equations  (5. 7) , (5. 9) , (5. 11) ,  and  (5.12) 
that  tracking  will  occur  in  the  sense  of  equation  (5.5) 
provided  only  that 

Zs  U  Zf  CV~  (5.31) 
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where  V~  is  the  open  unit  disc.  In  view  of  equations  (5.25), 
(5. 26)and(5.27) ,the  'slow'  and  'fast'  modes  will  satisfy 
the  tracking  requirement (5. 31) for  sufficiently  small  sampling 
periods  if  the  controller  and  transducer  matrices  Kq,  k^, 
and  M  are  chosen  such  that  both  Z^C^P  and  Z2  d  P  for  suf¬ 
ficiently  small  sampling  periods  and  Z^Op-  in  the  case  of 
plants  for  which  the  transducer  matrix  M  can  be  simultaneously 
chosen  so  as  to  satisfy  the  measurement  condition  expressed 
by  equation (5. 10) .  Moreover,  if  KQ  and  M  are  chosen  so  that 

■'*  A 

both  r (A)  and  T(A)  are  diagonal  transfer  function  matrices 
by  requiring  that 

F2B2Kq  =  (C2+MA12)B2K0  =  diag{a1,a2, . . . ,a^}  (5.32) 

where  1-a ^  €1  R  O  P  (j*l,2, . . . ,4)  in  the  case  of  T(A),  it 
follows  from  equation  (5. 24) that  the  transfer  function  matrix 
G( A)  of  the  discrete-time  tracking  system  will  assume  the 
diagonal  asymptotic  form  T ( A )  and  therefore  that  increasingly 
non-interacting  tracking  behaviour  will  occur  as  f  ->■ 
Furthermore,  such  tracking  behaviour  will  exhibit  high  accuracy 
in  the  face  of  plant-parameter  variations  provided  that  the 
steady-state  conditions  expressed  by  equation (5. 9) correspond 
to  'kinematic'  relations2iips  which  hold  between  the  state 
variables  as  a  consequence  of  the  fundamental  dynamical 
structure  of  the  plant. 

However,  increasingly  'tight'  tracking  behaviour  will 

not  in  general  occur  as  f  >  «  in  view  of  the  possible  presence 

in  the  plant  output  vector  of  'slow'  modes  corresponding  to 

the  poles  Z2  of  the  'slow'  transfer  function  matrix  r ( A) . 

But  the  presence  of  any  such  'slow'  modes  is  the  inevitable 
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consequence  of  the  introduction  of  appropriate  transducers 
which  generate  extra  plant  output  measurements  so  as  to 
ensure  that  rank  F2B2  *  ^  311(1  thus  render  plants  for  which 
rank  C2B2  <  i  amenable  to  fast-sampling  error-actuated  control. 
Indeed,  increasingly  'tight'  tracking  behaviour  is  achievable 
as  f  +  •  only  in  the  case  of  plants  for  which  rank  C2B2  =  £ 
(Bradshaw  and  Porter  19 8Q)  and  which  are  accordingly  amenable 
to  fast-sampling  error-actuated  control  without  the  necessity 
for  the  generation  of  extra  plant  output  measurements. 

5.4  Illustrative  Example 

These  general  results  can  be  conveniently  illustrated 
by  designing  a  fast-sampling  error-actuated  digital  flight 
controller  for  the  longitudinal  dynamics  of  an  aircraft 
governed  on  T  by  the  respective  state  and  output  equations 
(Kouvaritakis,  Murray,  and  MacFarlane  1979) 

"x^tfl  To  ,  1.132  ,  0  ,  0  ,  1 

x2(t)  0,0  ,  0  ,  1  ,  0  x2(t) 

x3(t)  =  0  ,  -0.1712  ,  -0.0538  ,  0  ,  0.0705  x3(t) 

x4(t)  0,0  ,  0.0485  ,  -0.8536  ,  -1.013  x4(t) 

x5(t)J  [p  •  0  >  -0.2909  ,  1.0532  ,  -0. 6859_J  [x5  ( t) 

0  ,  0  ,  o' 

0  ,  0  ,  0  fui(tf 

+  -0.0012  ,  1  ,  0  u2(t)  (5.33) 

0.4419  ,  0  ,  -1.6646  [u3 (t)_ 

0.1575  ,  0  ,  -0.0732 
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and 


V^t)' 

1,0,  0,0,0 

y2(t) 

m 

0  ,  0  ,  1  ,  0  ,  0 

0#l/0#0f0^ 

x^t) 
x2  (t) 
x3  (t) 
x4(t) 

X5  ( t)_ 


(5.34) 


from  which  it  can  be  readily  verified  that  the  first  Markov 
parameter 


0  ,  0  ,  ol 


C2B2  " 


-0.0012 

0 


1  ,  0| 

0  ,  oj 


(5.35) 


is  rank  defective.  In  case  {o^,a2,a^}  -  (1,1,1},  »  Kq, 

and 


0.25  ,  0 


M 


0  ,  0 
0  ,  0.25 


(5.36) 


it  follows  from  equations  (5.3),  (5.4),  and  (5.32)  that  the 
corresponding  transducers  and  fast-sampling  error-actuated 
digital  controller  are  governed  on  7  and  by  the  respective 
measurement  and  control- law  equations 


Wj^t)' 

“l  ,  0.283 

,  0  ,  0 

,  0.25' 

w2(t) 

as 

0 

0 

,  1  ,  0 

,  0 

w3(t)_ 

p  ,  1 

,  0  ,  0.25 

,  0  _ 

xx  ( t) 
x2(t) 
x3(t) 
x4(t) 


(5.37) 


lx5(t) 
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and 


u^kT) 
u2(kT) 
u3 (kT) 


=  f. 


-28.97  ,  0  ,  -1.274 

-0.0348  ,  1  ,  -0.001! 
-7.690  ,  0  ,  -2.740 

-28.97  ,  0  ,  -1.274 

-0.0348  ,  1  ,  -0.0015 


e1(kT)“ 

e2(kT) 

_e  3  ( kT )  _ 

■a 

z1(kT)' 

z2(kT) 

- 

,z3  (kT)_ 

V  (5.38) 


and  it  is  evident  from  equations  (5. 25) , (5. 26) , and (5. 27)  that 
Z]_  =  { 1-T ,  1-T ,  1-T  }  ,  Z2  -  {  1-4T ,  1-4T  ,  1-4T  }  ,  and  Z3  =  {0,0, 0>. 
It  is  also  evident  from  equation  (5.30)  that  the  asymptotic 
transfer  function  matrix  assumes  the  diagonal  form 


TU) 


4T 


A-1+4T  ' 


0 

0 


1 

'  I  ' 


0 

4T 


A-1+4T 


(5.39) 


and  therefore  that  the  closed-loop  discrete-time  flight-control 
system  for  the  longitudinal  dynamics  of  the  aircraft  will 
exhibit  increasingly  non-interacting  tracking  behaviour  as 
f  •  when  the  piecewise-constant  control  input  vector 
[ux(t)  ,  u2(t)  ,  u3(t)JT  =  [u^kT)  ,  u2(kT)  ,  u3(kT)]T, 
t  ^.[kT  ,  (k+l)T)  ,  kT  is  generated  by  the  fast-sampling 

digital  controller  governed  on  T T  by  equation  (5.38).  However, 
it  is  apparent  from  equation  (5.  3  5)  that  increasingly  ’tight’ 
tracking  behaviour  will  be  achieved  only  in  the  case  of  y2(t) 
in  view  of  the  presence  in  y1(t)  and  y3(t)  of  the  ’slow’ 
inodes  corresponding  to  the  poles  Z2. 
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The  actual  frequency-response  loci  G(eia)T)  for  tuT  jo  ,  2tt] 
are  shown  in  Figs  5. 1,  5. 2 ,  and  5. 3  when  T  «  0. 10,  0.05,  and  0.02, 
respectively,  and  it  is  clear  that  the  actual  frequency-response 
loci  approach  the  asymptotic  frequency-response  loci  r(eiajT) 
as  the  sampling  frequency  f  is  increased.  The  corresponding 
step-response  characteristics  ars  shown  in  Fig  5. 4 ,  Fig  5. 5 ,  and 
Fig  5.6  when  [v1(t)  ,  v2(t)  ,  v3(t)]T  =  [l  ,  0  ,  o]T,  [v^It)  , 
v2 (t)  /  v3(t)]T  «  [o  ,  1  ,  o]T,  and  [v^t)  ,  v2(t)  ,  v3(t)]T  = 

[o  ,  0  ,  l]T,  respectively,  and  it  is  evident  that  increasingly 
non-interacting  tracking  occurs  as  the  sampling  frequency  f 
is  increased  but  that  increasingly  'tight'  tracking  occurs 
only  in  the  case  of  y 2 ( t) . 

5,5  Conclusion 

Singular  perturbation  methods  have  been  used  to  exhibit 
the  asymptotic  structure  of  the  transfer  function  matrices 
of  discrete-time  tracking  systems  incorporating  linear 
multivariable  plants  which  are  amenable  to  fast-sampling 
error-actuated  digital  control  only  if  extra  plant  output 
measurements  are  generated  by  the  introduction  of  appropriate 
transducers  and  processed  by  inner-loop  compensators.  It 
has  been  shown  that  these  results  greatly  facilitate  the 
determination  of  controller  and  transducer  matrices  which 
ensure  that  the  closed-loop  behaviour  of  such  discrete-time 
tracking  systems  becomes  increasingly  non-interacting  as  the 
sampling  frequency  f  is  Increased.  These  general  results 
have  been  illustrated  by  the  presentation  of  the  frequency- 
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response  and  step-response  characteristics  of  a  discrete- 
time  flight-control  system  for  the  longitudinal  dynamics  of 
an  aircraft. 
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Fig  5 .1  (i)  :  G33(eia)T);  T  ®  0.10 
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Fig  5.3(a):  Gu  (elu,T)  ;  T  =  0.02 
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Fig  5 .4  ('>)  :  T  -  °  -  1° 
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Fig  5.6(h):  T  -  0.05 
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CONCLUSIONS  AND  RECOMMENDATIONS 
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6,1  Conclusions 


The  conceptual  and  computational  simplicity  of  singular 
perturbation  methods  (Porter  and  Shenton  1975)  in  the  design 
of  both  analogue  and  digital  tracking  systems  incorporating 
error-actuated  controllers  has  been  demonstrated.  This 
demonstration  has  been  effected  by  designing  such  controllers 
for  plants  which  have  previously  been  considered  by  other 
methods  (MacFarlane  and  Kouvaritakis  1977,  Kouvaritakis , 
Murray,  and  MacFarlane  1979)  so  that  detailed  comparisons 
between  alternative  design  techniques  can  readily  be  made, 
at  least  in  the  case  of  error-actuated  analogue  controllers. 

Moreover,  further  demonstration  of  the  power  of  the 
design  techniques  described  in  this  report  has  been  provided 
by  the  design  of  direct  digital  flight-mode  controllers  for 
the  F- 16  aircraft  (Porter  and  Bradshaw  1981).  In  addition, 
the  fact  that  high-gain  and  fast-sampling  error-actuated 
controllers  remain  highly  effective  in  the  presence  of  gross 
actuator  nonlinearities  has  been  conclusively  demonstrated 
(Porter  1981a,  Porter  1981b,  Porter  1981c,  Porter  1981d) . 
Thus,  the  same  basic  design  techniques  have  been  shown  to  be 
applicable  to  both  linear  and  nonlinear  multivariable  plants. 

6,2  Recommendations 

The  controllers  for  the  F-16  aircraft  designed  by 
singular  perturbation  methods  (Porter  and  Bradshaw  1981)  are 
extremely  robust  in  the  sense  that  the  same  controller 
remains  effective  throughout  a  wide  range  of  different  flight 
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conditions.  However,  in  case  even  higher  performance 
characteristics  are  required,  the  use  of  adaptive  controllers 
becomes  unavoidable  and  the  study  of  f ast-sampling  adaptive 
digital  controllers  is  therefore  recommended. 

Finally,  although  fast-sampling  digital  controllers  for 
the  F-16  aircraft  have  been  successfully  implemented  in  the 
Engineering  Dynamics  and  Control  Laboratory  at  the  University 
of  Salford  using  microprocessor  hardware  (Garis  1981) ,  the 
implementation  of  such  controllers  in  actual  aircraft  and  sub¬ 
sequent  flight  testing  is  recommended.  Such  practical  im¬ 
plementation  would  obviously  entail  close  liaison  with  the 
aerospace  industry,  based  initially  upon  the  wider  dissemination 
of  the  design  methodologies  described  in  this  report. 
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